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Abstract 

In this paper we calculate Kolmogorov's constant for magnetohy- 
drodynamic turbulence to one loop order in perturbation theory using 
the direct interaction approximation technique of Kraichnan. We have 
computed the constants for various E u (k)/E b (k), i.e., fluid to mag- 
netic energy ratios when the normalized cross helicity is zero. We 
find that K increases from 1.47 to 4.12 as we go from fully fluid case 
{E b = 0) to a situation when E u /E b = 0.5, then it decreases to 3.55 
in a fully magnetic limit (E u = 0). When E u /E b = 1, we find that 
K = 3.43. 

PACS number: 47.65. +a, 47.25.cg, 52.35 
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Kolmogorov M hypothesized that the energy spectrum E(k) of fluid tur- 
bulence in the inertial range (i.e., length scale intermediate between the en- 
ergy feeding and energy damping ranges) is a power law with a spectral index 
of —5/3, i.e., 

E(k) = K Ko ii 2 / 3 k- 5 /\ (i) 

where Kk is an universal constant called Kolmogorov's constant, k is the 
wavenumber, and II is the nonlinear cascade of energy which is also equal to 
the dissipation rate of the fluid. 

The calculations based on Direct interaction approximations 0, [J , renor- 
malization group technique [|J, self-consistent mode coupling || etc. show 
that Kko ~ 1.5 in three dimensions (3D). However, in two dimensions 
(2D) Kraichnan [[|, Olla M, and Nandy and Bhattacharjee |7| show that 
Kko ~ 6.4 in the region where inverse cascade of energy occurs. Experi- 
ments and numerical simulations || yield constants approximately same as 
those predicted by the above calculations. 

Magnetohydrodynamic (MHD) turbulence is more complex than the fluid 
turbulence. In MHD there are velocity field u and magnetic field B = B + b, 
where B is the mean magnetic field and b is the fluctuation in the magnetic 
field. It is customary to use Elsasser variables z 1 * 1 = u ± b/ y/A-np that can be 
interpreted as Alfven waves having positive and negative velocity and mag- 
netic field correlations respectively. Here p is the density of the magnetofluid. 
In this paper we assume incompressible approximation which in most of the 
situations correspond to a constant density of the fluid. Matthaeus and 
Zhou, Zhou and Matthaeus, and Marsch || constructed a phenomenology 
when z 3> B Q /y/4np, where Bq is the mean magnetic field or the magnetic 
field of the largest eddies. In this phenomenology, the energy spectra E^ik) 
of z ± fluctuations are proportional to k~ b ^ 3 and 

E ± (k) = K ± (n ± ) 4 / 3 (n^)- 2 / 3 A;- 5 / 3 , (2) 

where K ± are constants, which we will refer to as Kolmogorov's constants 
for MHD turbulence, and n ± (fc) are the nonlinear cascade rates of z fluctu- 
ations respectively. The normalized cross helicity a c , (strictly speaking spec- 
tral normalized cross helicity) defined as (E + (k) — E~ (k)) j \E + (k) + E~(k)), 
play an important role in MHD turbulence. In this paper, however, we re- 
strict ourselves to zero normalized cross helicity where, by symmetry, n + = 
n~ and K + = K~ = K. This choice of a c is motivated by the simplicity of 
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the calculation and by the fact that in the solar wind most of the fluctuations 
in the outer heliosphere (beyond earth) and some of them in the inner he- 
liosphere (from sun to earth) have negligible a c . There is another important 
energy called residual energy E R defined as the difference between the kinetic 
(E u ) and magnetic (E b ) energies. Note, however, that the residual energy is 
not conserved while both E + and E~ are conserved. 

In the other limit when z <C B /^4irp, Kraichnan and Dobrowolny et 



al. [10 1 argued that 



Tl + = Tl- = —E+(k)E-(k)k\ (3) 

where A is a constant called Kraichnan's constant. 

The solar wind observations indicate that the it exhibits MHD turbulence. 
The energy spectra of the solar wind has been found to be closer to fc _5//3 than 
to &;~ 3 / 2 [|Tlf| . However, it should be noted that these exponents are difficult 



to distinguish. Recently we have performed direct numerical simulations 
of MHD turbulence [0 in that we could not conclude whether the spectral 
index was 5/3 or 3/2. However, the nonlinear energy cascade rates IP 1 appear 
to follow Eq. (HP rather than Eq. (|3|). Motivated by the above mentioned solar 
wind observations, in this paper we assume that the energy spectra E ± {k) 
obey Kolmogorov-like phenomenology [Eq. (0)] for z ± > Bq. 

We determine for the first time the constant K associated with the Kol- 
mogorov spectrum in MHD turbulence. To this end we attempt to find K of 
Eq. (fj) using perturbative technique similar to DIA (see Fournier et al. and 
Camargo and Tasso |13| for application of Renormalization Group techniques 



to MHD). In this paper we only discuss the situations when E + {k) = E (k) 
and -Bo = 0. However, we vary E u /E b from to oo, i.e., from magnetic case 
to fully fluid case. The computation of when E~(k)/E + (k) is different 
from unity and when Bq is nonzero is under progress; forthcoming results 
will be presented elsewhere. 

The MHD equations in terms of z ± , in absence of B Q , are given by [|10J 

iztiKt) + v + k 2 zf^t) + V_k 2 Z?{\L,t) = ft 
+eM ijm (k) J dpzf (p, t)z±(k - p, t) 



where 



kjk m v i (? /(Attct^ 

Mij m (is.) = fcj\Pj m (k); -Pj m (k) = 5i m ; v± = - , (5) 
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v is kinematic viscosity, c is the speed of light, a is conductivity, ^ are 
forcing, and e is the expansion parameter which is set to one at the end. We 
assume independent forcing of z ± , i.e., df Sl /df S2 = 5 SlS2 , where Si,s 2 = ±. 
We solve for self-energy X (a 2 x 2 matrix), which is defined as 

G = G - G tG, (6) 



to first order using the technique of Wyld [|14] and Leslie || (the details of 
the calculation are be presented in ||15|| ). In our calculation we postulate the 
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frequency dependence of the z — z correlations functions as ||: 

where rj^ and r/^T are the reciprocal of the mean response time of z£ and z^ 
fluctuations respectively. The self-energy can be written as a following 2x2 
matrix: 

/ V"H — \ 1 / „4- „— \ 

(8) 



rf^ik) = lim ^k 2 J J dpdu' x 

[hfcp, q)(G^(p, J) + G^(p, c/))C*(k -p,cu- u') 
+b 2 (k, p, g)(G ±=F (p, u') + G^ip, u'))C*{k -p,u- a/) 
+6 3 (A;, p, g)(G ±± (p, c/) + G^(p, c/))C*(k - P, u - a/) 
+h(k, p, q)(G^(p, to') + G^F(p, ^))C i? (k - P, w - a;')] (9) 

where G is the Green's function, bi(k,p,q) = k~ 2 Bi(k,p,q), and B^s are 
the functions of (x, y, z), the cosines of the angles between (p,q), (q,k), and 
(k,p) respectively; B^s given in the appendix of Leslie ||. Substitution of X 
of (H) in Eq. (||) and ignoring kinematic viscosity and resistivity as compared 
to S/A; 2 yields G 



G{k,uj) 



1 f + a 



2 2 



Note that when C + (fc) = C (A;) = C R (k), i.e., when the magnetic energy is 
zero, the above equations yield r] + = rf~ = rjfi ui d- 
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As mentioned earlier, in this paper we restrict ourselves to situations when 
C + (k) = C~(k) = C(k) and assume the Kolmogorov scaling law, i.e., E(k) 
are proportional to k~ 5 ^ 3 . In this situation rf^ = rj^ by symmetry. Since the 
energy spectra E(k), which is equal to C (k) / (Auk 2 ) , is proportional to fc~ 5//3 , 
the correlation functions C [k) and self-energies will be || 

C (k) = l-KYi 2 ' 3 ^ 11 ' 3 ; jj^k) = AU^ 3 k 2 / 3 , (10) 

where A is a constant. We denote the ratio C R /C by a. Using the above 
quantities and with the change of variables p = (k and q = nk we obtain 

x [a(6i(fc,p,g) + b 4 (k,p,q)) + b 2 (k,p,q) + b 3 (k,p,q)] . 



The integrals of Eq. fll2|) suffer from the well known "infrared problem" which 
comes from the strong dynamic coupling of fluctuations with widely differ- 
ing wavenumbers. Over the years, various techniques have been developed 
to tackle the difficulty in the context of pure fluid turbulence. The ear- 
lier methods (essentially introduction of a cut off) are discussed in LeslieQ. 
Later work involving either Renormalization group technique [Q, or the La- 
grangian or semi-Lagrangian pictures ||16|| , or self-consistent screening || 



shows how the theory can be made naturally finite and cut off independent. 
The infrarred difficulties associated with Eq. ( |T2"D can be similarly avoided. 
Knowing that the full theory is constrained to be finite, we adopt the practi- 
cal procedure of evaluating the integral in Eq. ( |I2"D with a cut off k = Xk and 
choosing A = 1 so that the pure fluid value of A 2 / K are obtained correctly 
when a = 1. Thereafter A is not varied. With this procedure we obtain 
A 2 /K as a function of a. For some of the characteristic values of a, A 2 /K is 
listed in Table 1. 

To obtain the numerical value of K we need another equation involving 
A and K. To this end, we derive an expression for energy cascade rate 
n in terms of correlations C(k) and Green's functions G(k). Following an 
approach similar to that of Leslie 0, when C + {k) = C~(k) = C(k), we 
obtain 

(^- t + 2u + k 2 ^j E(k) + 2v_k 2 E R (k) = T(k, t). (12) 
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Using T(k, t) we can write an expression for energy cascade rate, which is 
II = — J dk'T(k',t). Substitution of C's and 77's, and change of variables 
k' = k/u,p = v(k/u), and q = w(k/u) yield an expression for A/K 2 , that is, 



A r 1 1 f 1+v 

— = dvln- dw[y(l,v,w) + ty(l,w,v)] (13) 
K Jo v Jv* 



where 



V(l,v,w) = ^vw(l + v 2/s + w 2/3 )~ 1 (14) 

x [6x(l, v, w)w- u/3 ((l + a>- 11/3 - 2a) 

+ (1 + a 2 )b 2 (l, w, v)v- u/3 (w- 11/3 - 1) 

- 2ab 3 (l, v, w)w- u/3 - (1 + a 2 )h(l, v, w)v' u/3 ], 

and v* = max(f , |1 — v\). See Table 1 for A/K 2 for various values of a. 

We calculate K from A 2 / K and A/K 2 . The values of K for various values 
of a are listed in Table 1. We have also calculated the constant K in 2D. 
In 2D when E u = E b , we choose A = 1, same as A of 3D, but not 0.065 
which yields K = 6.6 in fluid limit. The choice of same A in 2D MHD was 
motivated by fact that in MHD turbulence, forward cascade of energy and 
inverse cascade of cross helicity occur in both 2D and 3D|]T7|. Note that the 



behaviour of fluids turbulence in 2D and in 3D are dramatically different. 

The results presented in this paper are in good agreement with the 
simulation results of Verma et al. []T2] (K = 3.7 in three-dimensions for 
E~/E + = 0.6, and average K = 6.6 in two-dimensions for E~ / E + = 1,0.6). 
In 3D simulation, however, only a single run with a relatively lower resolution 
(128 3 ) was performed. We need to perform more runs to come to a definite 
conclusion. 

We find that K monotonically increases from 1.47 to 4.12 as we go from 
fully fluid case (E R = E) to a situation when E R /E = -1/3, i.e., E u /E b = 
0.5, then it decreases and finally reaches 3.55 in fully magnetic limit (E R = 
—E). The maxima of K occurring for E u / E b = 0.5 is a curious result because 
average E u /E b for the solar wind in the outer heliosphere, where E + m E~ 
is a good approximation, is 0.5. Note that maximum K corresponds to 
minimum energy dissipation, hence it appears as if the solar wind is settling 
to a minimum dissipation state. Further investigation in this direction could 
possibly yield interesting results regarding minimum dissipation states. 
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Recent temperature evolution studies of solar wind plasma show that 
turbulent heating in the solar wind can account for the observed temperature 
evolution when K = 1.0 for nonAlfvenic streams (streams for which E + pa 
E~) In the inner heliosphere, where the solar wind is fluid dominated, 
the turbulent heating is possibly the major heating mechanism. However, in 
the outer heliosphere, since K is in the range of 3 to 4, turbulent heating could 
account for only 20-25 % of the total dissipation; rest of the dissipation should 
be provided by other sources, e.g., shock heating, stream-stream interactions 
etc. 

To conclude, we have calculated the Kolmogorov's constant for MHD 
turbulence for zero a c to one loop order using perturbative technique similar 
to the direct interaction approximation of Kraichnan. We find that in fluid 
dominated case, K pa 1.5—2, but as magnetic energy increases, K increases as 
well and reaches a maximum value of 4.12 when E b = 2E U , then it decreases 
and reaches 3.55 for a fully magnetic case. These results are consistent with 
the recently performed simulation results and also shed light on the evolution 
of E u /E b ratio and temperature in the solar wind. 

We thank Malay Nandy for fruitful discussion and R. Sethuraman for 
help in computation. We also thank one of the referees for careful reading of 
the manuscript and for useful comments. 



7 



References 

[1] A. N. Kolmogorov, C. R. Acad. Sci. URSS, 30, 301, 1941. 
[2] R. H. Kraichnan, J. Fluid Mech., 47, 525, 1971. 

[3] D. C. Leslie, Developments in the theory of Turbulence, Claredon, Ox- 
ford, 1973. 

[4] V. Yakhot and S. A. Orszag, J. Sci. Comput., 1, 3, 1986; D. Ronis, 
Phys. Rev. A, 36, 3322, 1987; D. Forster, D. Nelson, and M. Stephen, 
Phys. Rev. A, 16, 732, 1977; For recent calculations and reviews refer to 
articles in Transition, Turbulence, and Combustion, ed. M. Y. Hussaini 
et al., Kluwer, 1994. 

[5] J. K. Bhattacharjee, Phys. Fluids, A 3, 879, 1991; V. I. Belinicher and 
V. S. L'vov, Sov. Phys. JETP, 66, 303, 1987. 

[6] P. Olla, Phys. Rev. Lett., 67, 2465, 1991. 

[7] M. K. Nandy and J. K. Bhattacharjee, in press, Int. J. Mod. Phys., 
1995. 

[8] H. L. Grant, R. W. Stewart, and A. Molliet, J. fluid Mech., 12, 241, 
1962; U. Frisch and P. Sulem, Phys. Fluids, 27, 1921, 1984. 

[9] W. H. Matthaeus and Y. Zhou, Phys. Fluids, Bl, 1929, 1989; Y. Zhou 
and W. H. Matthaeus, J. Geophys. Res., 95, 14881, 1990; E. Marsch, 
in Reviews of Modern Astronomy, Vol. 4, ed. G. Klare, Springer- Verlog, 
1990. 

[10] R. H. Kraichnan, Phys. Fluids, 8, 1385, 1965; M. Dobrowolny, A. Man- 
geney, and P. Veltri, Phys. Rev. Lett., 45, 144, (1980). 

[11] W. H. Matthaeus and M. L. Goldstein, J. Geophys. Res., 87, 6011, 1982; 
E. Marsch and C.-Y. Tu, J. Geophys. Res., 95, 8211, 1990. 

[12] M. K. Verma, Ph. D. thesis, University of Maryland, College Park, 1994; 
M. K. Verma, D. A. Roberts, M. L. Goldstein, S. Ghosh, and W. T. 
Stribling, preprint, 1995. 



8 



[13] J. D. Fourrnier, P. L. Sulem, and A. Pouquet, J. Phys. A, 15, 1393, 
1982; S. J. Camargo and H. Tasso, Phys. Fluids B, 4(5), 1199, 1992. 

[14] H. W. Wyld, Ann. Phys. (NY), 14, 143, 1961. 

[15] M. K. Verma and J. K. Bhattacharjee, to be sent to Phys. Rev. E, 1995 
[16] R. H. Kraichnan, Phys. Fluids, 8, 575, 1965. 

[17] A. C. Ting, W. H. Matthaeus, and D. Montgomery, Phys. Fluids, 29, 
9695, 1986. 

[18] M. K. Verma, in press, J. Geophys. Res., 1995. 



9 



E H /E 


E u /E b 


A 2 /K 


A/K' 2 


K 


1 


oo 


0.114 


0.190 


1.47 





1 


0.176 


0.066 


3.43 


-1/3 


1/2 


0.197 


0.053 


4.12 


-2/3 


1/5 


0.218 


0.056 


4.11 


-1 





0.239 


0.073 


3.55 


Q 2V 


1 


0.070 


0.021 


5.41 



Table 1: Kolmogorov's constant K of MHD turbulence for various values of 
a = E R /E (assuming E + = E~). The symbols E,E R ,E u ,E b denote the 
total energy, residual energy, kinetic and magnetic energy respectively. 
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